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ABSTRACT 

We report high angular resolution observations of the HCN (3-2) line emission in the circumstellar 
envelope of the 0-rich star W Hya with the Submillimeter Array. The proximity of this star allows us 
to image its molecular envelope with a spatial resolution of just ~ 40 AU, corresponding to about 10 
times the stellar diameter. We resolve the HCN (3-2) emission and find that it is centrally peaked and 
has a roughly spherically symmetrical distribution. This shows that HCN is formed in the innermost 
region of the envelope (within ~ 10 stellar radii), which is consistent with predictions from pulsation- 
driven shock chemistry models, and rules out the scenario in which HCN forms through photochemical 
reactions in the outer envelope. Our model suggests that the envelope decreases steeply in temperature 
and increases smoothly in velocity with radius, inconsistent with the standard model for mass-loss 
driven by radiative pressure on dust grains. We detect a velocity gradient of ~ 5 km s _1 in the 
NW-SE direction over the central 40 AU. This velocity gradient is reminescent of that seen in OH 
maser lines, and could be caused by the rotation of the envelope or by a weak bipolar outflow. 
Subject headings: circumstellar matter — radio lines: stars — stars: individual (W Hya) — stars: 
late-type — stars: winds, outflows 



1. INTRODUCTION 

The detection at radio frequencies of various molecules 
in the circumstellar envelopes of evolved stars has con- 
siderably improved our understanding of circumstellar 
chemistry. This is particularly true for the prototype 
of carbon stars, IRC+10216, on which much atten- 
tion has been focused. Several dozens of molecules 
have been detected in its circumstellar envelope (e.g., 
ICernicharo et al.ll2000l ). reflecting its complex and rich 
chemistry. The inventory of molecules together with 
high angular resolution images of their spatial distribu- 
tion have led to the construction of sophisticated chemi- 
cal m odels for the circumstellar envelope of IRC+ 10216 
(e.g., iGlassgold et"aII[T987l : iMillar fc Herbstl [l99l that 
now form the basis of our understanding of the chem- 
istry in C-rich envelopes. 

The chemistry in O-rich envelopes, on the other 
hand, was not expected to be as rich because equi- 
librium chemistry calculations indicate that nearly all 
the carb on should b e locked into CO and nitrogen 
into N 2 ([Tsujil Il973h . Indeed, before 1985, the in- 
ventory of molecules detected around oxygen stars was 
sparse and di d not include carbon-b e aring molecules, ex- 
cept of CO. iDeguchi fc Goldsmith! ((1985) reported the 
first detection of HCN towards O-rich sta rs. The list 
of HCN detections was later extended by I Jewell et al.l 
(198 61). iLindayist et all (119881). iNercessian et all (ll989f). 
l01ofssonet_alJ (I1998D an d iBieging et all (|2000h . Re- 
cently, IZiurvs et al. (2007) also emphasized the chemical 
complexity in the envelope of the O-rich supergiant star 
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VY CMa, following the detection of various molecules 
such as HCO+, CS, NaCl, PN and SiS. 

Two competing models have been proposed to account 
for the formation of HCN in O-rich en yelopes: photo- 
chem ical reactions in the outer envelope (|Charnley et al.l 
Il995[ ) , or gas-phase non-equilibrium chemical reactions in 
the inner region close to the stellar photosphere (within 
few tens of AU fro m the star) due to shocks driven 
by st ellar pulsation (|Duari et al.lll999l : iDuari fc Hatchelll 
2000). The two models predict very different HCN spa- 
tial distribution: in a hollow shell-like structure of radius 
200 to 1000 AU depending on the mass loss for photo- 
chemical reactions, or a centrally concentrated distribu- 
tion extending up to the photodissociation radius of HCN 
for the shock-driven chemical reactions. 

Current observational data seem to favor pulsation- 
driven shock chemist ry formation of HCN in O-rich stars. 
Bieg ing et al.l (|2000l ) conducted single-dish observations 
toward a sample of 16 O-rich stars in the HCN (3-2) and 
(4-3) lines. The detection of these high density tracers 
is inconsistent with photochemical production of HCN 
in the outer envelope. The detection of the HCN (3-2) 
(0,l lc ,0) and (8-7) transitions toward \ Cyg gave further 
support for the formation of HCN in the innermost re- 
gion of the envelope, within ~ 20 stell ar radii (i.e. ~ 30 
AU), b ased on excita tion arguments (IDuari fc Hatched 
120001 ). iMarvell ([2005) reported the first interferometric 
observations of HCN in the circumstellar envelopes of 
the O-rich stars IK Tau and TX Cam. In both stars, the 
HCN emission appears to be largely concentrated within 
^ 800 AU, as predicted by shock chemistry models. 

If HCN indeed forms close to the stellar photosphere, 
it should be an excellent tracer of the inner envelope, and 
perhaps even of the wind acceleration zone where, in the 
standard model, dust particles form and radiation pres- 
sure on dust becomes most effective. Maser (e.g., H2O 
or OH) emission in the inner envelope can be observed 
with very long baseline interferometry, but traces only 



2 



Muller et al. 2008 



a limited radial range. On the other hand, there are so 
far only few reports of high angular resolution observa- 
tions of other therma l lines in the inner envelope, such 
as Si O (2-1) v = dLucas et all 119921: ISahai fe Biegind 
1993; Schoier et al. 2004), none of them with angular res- 
olution higher than ~ 100 AU. 

In this letter, we report very high spatial resolution 
(~ 40 AU) observations of the HCN (3-2) emission to- 
ward W Hya with the Submillimeter Array. W Hya is 
one of the closest O-rich stars at a distance o f 78 pc 
(b ased on Hipparcos d ata JPerrvman et al]ll997l revised 
bv lKnapp et ^aLll2QQ3[ ). and shows the brightest HCN (3- 
2) and (4-3) emission amongst th e sample of O-rich stars 
observed by Biegi ng et al.l (|2000l ). It is therefore an excel- 
lent target in which to investigate the distribution and 
kinematics of HCN in the circumstellar envelope of an 
O-rich star. 

2. OBSERVATIONS 

We observed W Hya on 2008 April 13 and 15 with the 
Submillimeter Array 4 (SMA) in its very extended con- 
figuration. All eight antennas were operating on the first 
day of our observations, and six on the second day. The 
zenith atmospheric opacity was ~ 0.2 at 225 GHz dur- 
ing both nights, resulting in system temperatures that 
changed between 200 and 600 K depending on source 
elevation. 

The heterodyne SIS receivers were tuned to the fre- 
quency of the HCN (3-2) transition at 265.886 GHz in the 
lower side band. The correlator was configured to give a 
spectral resolution of 0.8125 MHz, which corresponds to 
a velocity resolution of ~ 0.9 km s _1 . The line-free chan- 
nels in both the upper and lower side bands, spanning a 
total bandwidth of ~ 3.8 GHz at 1.1 mm, were used 
to make a map of the continuum. The bandpass of the 
individual antennas was derived from the bright quasar 
3C273. Flux calibration was derived by observing Ti- 
tan and Callisto. The quasars 1334—127 and 1313—333, 
located within 16° of W Hya, were observed every 15 
min for complex gain calibration. Data reduction was 
done separately for both tracks using MIR/IDL. The cal- 
ibrated visibilities were written out in FITS format, and 
then converted into GILDAS format for imaging purpose 
using MAPPING. 

The projected baselines ranged from 22 to 510 m. Nat- 
ural weighting of the visibilities yielded a synthesized 
beam of 0.55" x 0.40" at a position angle of 26°. We used 
the CLARK algorithm to deconvolve the image. For the 
HCN (3-2) channel map, we estimate the la rms noise 
level to be 15 mJy/beam in each 0.9 km s _1 channels. For 
the 1.1 mm continuum, the noise level is 3 mJy/beam. 
The primary beam of the SMA antennas spans a FWHM 
of 46.5" at 266 GHz. 

3. RESULTS 

The HCN (3-2) channel maps are shown in Figure [1] 
together with the 1.1 mm continuum map. The molecu- 
lar emission, which is clearly resolved, is centrally peaked 
and has a roughly spherically symmetrical distribution. 
The continuum emission is unresolved and has a total 

4 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



flux density of 270 ± 15 mJy at 1.1 mm, consistent with 
the va lue of 280 ±30 mJy obtained by Ivan der Veen et al.l 
(1995) with the James Clerk Maxwell Telescope. We esti- 
mate a stellar blackbody contribution of 225 mJy based 
on the stellar parameters given in Table 1, and which 
therefore dominates the 1.1 mm continuum emission. 

W Hya is known to have a very extended dusty enve- 
lope with a radius of ~ 30' (Hawkins ll990t ). The pho- 
todissociation radius of HCN, however, is expected to be 
much smaller, of ord er of 1", based on computations by 
lOlofsson et al.l (|l998l ) (see their Eq.2). From a Gaussian 
fit of the visibilities, we find that the HCN (3-2) emitting 
region has a radius at FWH M of about -6", c ompara- 
ble to the size measured by iLucas et al.l (|1992D for the 
SiO (2-1) emission. Compa ring; our data with th at from 
single-dish observations by iBieging et al.l (|2000l ), we es- 
timate that we recover > 70% of the total line intensity 
within 28" of the star. A spectrum of the HCN (3-2) 
line, extracted at the peak position of the emission, is 
shown in Figure [2j The peak brightness temperature 
of the line is ~ 300 K. The line profile is clearly non- 
symmet ric with respect to the s ystemic velocity of ~ 40 
kms" 1 ([Cernicharo et al.lll997h . 

A visual inspection of the channel maps reveals a weak 
velocity gradient in the NW to SE direction. A fit of 
the visibilities yields a velocity gradient of ~ 5 km s _1 
at a position angle of 125° over the 0.5" central region. 
Position-velocity diagrams along this direction, as well 
as on a perpendicular axis, are shown in Figure [3j 

4. MODELING 

The peak brightness temperature, line profile and ra- 
dial distribution of the HCN (3-2) gas provide very useful 
constraints on the wind velocity and the physical con- 
ditions of the gas in the vicinity of W Hya. We used 
the molecular excitation and radiative transfer code of 
iDinh-V-Trung fe Nguyen-Q-Rieul (|2000f ) to interpret our 
observations and other relevant data. We assumed a 
spherical symmetry and did not attempt to reproduce 
the weak velocity gradient. The model parameters are 
summarized in Table 1 and reproduce remarkably well 
our data (see Figj2]&[4]). 

In the inner region of the envelope, the HCN molecules 
are excited by collisions with H2 and by the absorption 
of stellar photons at 3 fim (to ^3 = 1 vibrational state), 
7 jam (to 02°0 vibrational state), and 14 jam (to 01 1 
vibrational state). In the case of W Hya, where little ab- 
sorption by dust is expected, the pumping route through 
the z/3 = 1 and 01 1 states are the most important be- 
cause of the strong stellar radiation field at the corre- 
sponding wavelengths and large transition rates. In our 
code, we take into account the vibrational ground state, 
(01 1 0) and ^3 = 1 vibrational states, and include explic- 
itly all hyperfine levels up to J = 15. In W Hya, we 
find that the vs = 1 absorption band of HCN falls withi n 
the deep H2O absorption band (|Justtanont et al. 2004). 
The water molecules responsible for this absorption band 
are hot and expected to be very close to the central star, 
presumably interior to the HCN emitting region. Thus, 
for the sake of simplicity, we reduced the stellar radiation 
field by 50% at 3 fim. 

The mass loss rate o f W Hya is still quite uncertain. 
IZubko fc Elitzurl (|2000f ) derived a mass loss rate of about 
10 -6 M yr _1 (scaled to the revised distance of 78 pc), 
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while lJusttanont et al.l (|2QQ5l ) favored a much lower rate 
of 2-3 x 1(T 7 M yr- 1 . We adopt here a value of 5 x KT 7 
M yr -1 . The abundance of HCN relative to molecular 
hydr ogen predicted from non-equilibrium c h emica l mod- 
els (jDuari et all (|1999j ): iDuari fc Hatchelll (|2000h ) may 
vary between a few times 10 -6 and a few times 10 -5 , 
depending on the shock velocity. We adopt a represen- 
tative constant value of 10 -6 . We note, however, that 
this assumption leads to a very slight overestimate of 
the HCN (3-2) intensity in the outer part of the HCN 
envelope (see FigHJ. Most likely, the abundance of HCN 
is gradually decreasing with radius due to the photodis- 
sociation by external radiation field. We also assumed a 
local turbulent velocity of 1 km s _1 in the envelope. 

We tried various expansion velocity laws and found 
that the relation V(r) ex log(r), starting from 2 km s _1 
at an inner radius of ~ 10 14 cm and reaching 7 km s _1 
at a radius of 10 16 cm, produces a satisfactory fit to 
the HCN (3-2) line pr ofile, and is con si stent with OH 
maser observations of ISzymczak et al.l (|1998f ). Other 
faste r rising velocity laws, such a s those prescribed 
bvlDeguchi fc Nguyen-Q-Rieul (fl990h or lZubko fc Elitzurl 
(2000), do not reproduce the rounded-top and asymmet- 
ric profile of HCN (3-2) as seen by the SMA (they actu- 
ally produce double-peak profiles). 

The rounded-top line profile indicates that the HCN 
(3-2) is optically thick. In addition, we find that the 
excitation temperature of the HCN (3-2) transition fol- 
lows closely the kinetic temperature of the gas. Un- 
less the wind is highly clumpy, resulting in small fill- 
ing factor in the beam, the brightness temperature of 
the HCN (3-2) emission, resolved by the SMA, is there- 
fore directly related to the temperature profile in the 
inner region. The very high temp erature profile used by 
iDeguchi fc Nguyen -Q-Rieu (1990) or the one designated 
as GS in Figure 3 of IZubko fc Elitzurl (|2000f ) are clearly 
inconsistent with our data. The shallowe r temperature 
profile derived bv lZubko fc Elitzurl (|2000l ) also produces 
too high a brightness temperature for the HCN (3-2) line 
because the kinetic temperature is above 400 K within 
an inner region of 50 AU in diameter, which is compa- 
rable to the SMA synthesized beam of ~ 0.5". In our 
model, we need to use a steep temperature profile T(r) 
oc r _1 , starting with To = 650 K at the base of the en- 
velope. The predicted line profile is compared with that 
observed in Figure O Our model also provides a good 
match to the single-di sh profiles of the HC N (3-2) and 
(4-3) lines observed by Biegi nlTet al.l (|2000l ). 

We note that for the HCN (3-2) line, the resulting pro- 
file is asymmetric and does not peak at the systemic ve- 
locity. That is due to the strong self-absorption arising 
from the combination of steep temperature profile and 
slow acceleration in the inner envelope. This effect can 
also explain the difference in the inferred systemic ve- 
locity from Gaussian fitting of opti cally thick lines suc h 
as HCN (4-3), SiO(5-4) and (8-7) (|Bieging et all 120001 ). 
with respect to other lin es such as SiO(2-l) and CO lines 
(|Cernicharo et al1ll997l ). 

5. DISCUSSION AND PERSPECTIVES 

In the model proposed bv lCharnlev et all (|1995j ). HCN 
is formed in the outer envelope as a result of photochem- 
ical reactions initiated from methane. The distribution 
of HCN is predicted to peak at a radius > 200 AU for a 
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Fig. 1. — 1.1 mm continuum emission (top left corner box) and 
HCN (3-2) channel maps of W Hya as observed with the SMA very 
extended configuration. Contours are every 25 mJy/beam for the 
continuum map, and every 0.75 Jy/beam (5cr), corresponding to 
60 K for the synthesized beam of 0.55" x 0.40" (P.A.= 26°), for 
the channel maps. The white cross indicates the position of the 
phase center, set at the stellar position from the Hipparcos catalog 
in the year 2000 (R.A. J20 oo = 13 h 49 m 01!998 and Dec. J2 ooo = 
-28°22 / 03( / 49). The white dot gives the expected position of the 
star at the date of our observations in 2008, due to its proper 
motion. The small offset (< 0.15") between this position and the 
peak emission may be due to baseline errors. 
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Fig. 2. — Brightness temperature at the peak of the HCN (3-2) 
emission as a function of velocity, with model prediction overlaid 
(dashed line). Note that the blue-shifted part of the spectrum is 
strongly affected by self-absorption. 
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Fig. 3. — Posit ion- velocity diagrams along two perpendicular 
axes. Note the small velocity gradient along the SE-NW direction. 
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TABLE 1: MODEL PARAMETERS 



Distance 

Systemic velocity (LSR) 
Effective temperature 
Stellar radius 

Inner radius 
Outer radius 
Mass-loss rate 
Temperature profile 
Launching velocity 
Terminal velocity 
Velocity law (r < Rout) V(r) 
Abundance of HCN 
Local turbulent velocity 



STAR 
D = 78 pc 
Vsys = 40.4 km s _1 
T eff; * = 2500 K 
R* = 2.73 x 10 13 cm 
ENVELOPE 
R in = 1 x 10 14 cm 
R out = 3 x 10 15 cm 
M = 5 x 10- 7 M yr- 1 
T(r) = 650 K (r/R^)- 1 
V = 2 km s- 1 
Voo = 7 km s _1 

(Voc-Vo)loglo(r/R in ) 
log 10 (10 16 cm/R in ) 

[HCN]/[H 2 ] = 10- 6 
^turb = 1 km s" 1 



(1) 
(2) 
(3) 
(4) 



(5) 
(4) 



Vo 



Referen ces. — ( 1) Knap p et al 
(fl997h: (3) Haniffetal. (1995ft; (4) 
IMivoshi et al. (1994). 



d2003h: f 2) [C ernichar o et al.l 
Justtanont et al. (2005|); (5) 



Radius (arcsec.) 

Fig. 4. — Azimuthal average of the HCN (3-2) emission as a 
function of radius at Vlsr = 42.1 km s _1 , corresponding to the 
peak of the line profile in Figure [2] Our SMA data are indicated 
by open squares with error bars overlaid. The model prediction is 
shown in dashed line. 

mass-loss rate > 10 _T M yr -1 . Our SMA observations, 
however, clearly reveal that in the case of W Hya, HCN is 
present, and thus should form, much closer to the stellar 
photosphere, within 20 AU. This is consistent w i th the 
shoc k chemistry mode l deve loped by lDuari et al.l ([l999) 
and lDuari fc Hatchelll (|2QQQh . In this model, stellar pul- 
sations induce strong and periodic shocks in a narrow re- 
gion above the photosphere. Thermal equilibrium (TE) 
abundances are assumed as initial conditions in the pho- 
tosphere, and the evolution of chemical abundances is 
investigated in the post-shock region. As a net result of 
shock chemistry, several carbon-bearing species, such as 
HCN, CS and CO2, are produced in significant amounts. 
Noticeably, the abundance of HCN can be increased by 



about 5 orders of magnitude with respect to its former 
TE value. The route to HCN is through the reaction of 
CN with H2, which is sensitive to temperature and be- 
comes very efficient in the post-shock layer. Since HCN 
is chemically stable and does not participate in the for- 
mation of dust grains (unlike, e.g., SiO) in the inner en- 
velope, it further remains unaltered as it travels through 
the envelope, until it reaches the photodissociation region 
of the outer wind. As a result, HCN is a good tracer of 
the kinematics in the envelope. 

The requirement of a smoothly increasing velocity with 
radius to reproduce the HCN (3-2) line profile for W Hya 
and maser observations (|Szymczak et al.ll998j ) are incon- 
sistent with a fast acceleration of the wind within < 10 
stellar radii, as predicted by models of mas s- loss driven 
by radiative pressure on d ust grains (e.g., lKwokl fi975: 
iGoldreich fc Scovilld 1 1 9 76h . Observations in other lines 
tracing different radii also called for a slow wind accel- 
eration, which could be explained if dust grain forma- 
tion occirrs_over large extent up to several hundreds of 
AU (|Lucas et al. 1992) , and/or if grain propertie s change 

Cohenlll986h . Alter- 



(jijucas ex ai.m yyzK ana / or 11 \ 
through the envelope ([Chapman ■ 
natively, the slow acceleration could result from the in- 
effici ency of dust drag when the mass loss is small (see, 
e.g.. lBowenlll988UWinters et al.ll2000D . 

The velocity gradient along the SE-NW axis (see 
Fig El) is about 5 km s _1 over 0.5". MERLIN o bserva- 
tions of OH maser lines bv lSzvmczak et aTl (|l998l ) reveal 
a remarkably simil ar trend in ampl i tude a nd position an- 
gle. Interestingly, lLattanzi et al.l (|l997l ) measured the 
visible diameter of W Hya along two orthogonal direc- 
tions and found a significant asymmetry of 20%, with 
elongation in a direction comparable to that of the ve- 
locity gradient seen in the envelope. It is not clear, how- 
ever, whether this apparent elongation is due to some 
oblateness of the star or to a nonsymmetric brightness 
distribution over the stellar disk. The velocity gradient 
that we observe in the inner envelope of W Hya could be 
either due to the rotation of the envelope or to a weak 
bipolar outflow. 

The next jump in angular resolution achievable with 
the Atacama Large Millimeter Array will allow us to 
probe the envelope even deeper, possibly down to the 
dust formation zone. Such observations could directly 
test the validity of shock chemistry models. In particu- 
lar, very high angular resolution observations at differ- 
ent phases of the stellar variablity could directly show 
the effects of the pulsations on the formation of dust and 
mole cules, and possible modulat ions of the mass loss (see, 
e.g.. Diamond fc Kemba"ill2003l ). 
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